INT-PUB-11-055 



On the Role of Charmed Meson Loops in Charmonium Decays 

Thomas Mehen and Di-Lun Yang 1 

1 Department of Physics, Duke University, Durham, NC 27708, USA 

(Dated: November 22, 2011) 

Abstract 

We investigate the effect of intermediate charmed meson loops on the Ml radiative decays J /if) —> 
7/ c 7 and if)' — > i]P^f as well as the isospin violating hadronic decays if/ — > J /if) ir° '(n) using heavy 
hadron chiral perturbation theory (HHxPT). The calculations include tree level as well as one loop 
diagrams and are compared to the latest data from CLEO and BES-III. Our fit constrains the 
couplings of IS and 2S charmonium multiplets to charmed mesons, denoted g 2 and g' 2 , respectively. 
We find that there are two sets of solutions for g 2 and g 2 . One set, which agrees with previous 
values of the product 5252 extracted from analyses that consider only loop contributions to if/ —> 
J /if) 7T° (77) , can only fit data on radiative decays with fine-tuned cancellations between tree level 
diagrams and loops in that process. The other solution for g 2 and g 2 leads to couplings that are 
smaller by a factor of 2.3. In this case tree level and loop contributions are of comparable size 
and the numerical values of the tree level contributions to radiative decays are consistent with 
estimates based on the quark model as well as non-relativistic QCD (NRQCD). This result shows 
that tree level HHxPT couplings are as important as the one loop graphs with charmed mesons in 
these charmonium decays. The couplings g 2 and g 2 are also important for the calculations of the 
decays of charmed meson bound states, such as the X(3872), to conventional charmonia. 
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Many of the static properties and decays of charmonium states can be understood within 
a framework in which these states are viewed as non-relativist ic bound states of charm and 
anticharm quarks. This includes the quark model [l| as well as the modern QCD-based 
approach of non-relativistic QCD (NRQCD) |2j, which allows for systematic treatment of 
charmonium properties as an expansion in a s and v c , where v c is the relative velocity of the 
charm-anticharm quarks. Despite many successes there remain specific transitions that are 
not well understood quantitatively. Examples of decays that are not completely understood 
are the hadronic decays ip' — > J/ip(7T°,r]), and the radiative decays to J/ip — > 7/ c 7 and 
ip' — > rjc 7. The hadronic decays violate isospin, in the case of a final state with 7r°, or 
SU(3), when the final state is rj. As a consequence the ratio of this decay is sensitive to light 
quark masses ^, J]. The value of the light quark mass ratio extracted from the measured 
decay rates la], m u /md = 0.4 ± 0.01, differs significantly from the result extracted from 
meson masses in chiral perturbation theory, m u /m^ = 0.56 t[ S|- For the radiative decays 
the experimentally measured rates differ from quark model expectations. For example, a 
non-relativistic quark model calculation of J ftp — > rjcj (ip' — > ?7c7) yields a prediction of 
ft 3 (ft 0) keV, 1 whereas the experimental results are 1.57 ± 0.38 (0.97 ± 0.14) keV [10]. 



In Ref. 



111 ]. NRQCD is used to analyze the decay J/x/j — > rjc'-f, and the authors show that 



0(v 2 ) corrections can lower the rate so that the theoretical prediction is consistent with data. 



However, no attempt has been made to understand radiative decays of ip' in this 
For reviews of these puzzles and others in charmonium physics, see Refs. {9], I12I 



ramework. 



13|. 



Recently, Ref. 14| proposed that the hadronic decays mentioned above are dominated 



by loop diagrams with virtual D mesons 



he decays are calculated using Heavy Hadron 



Chiral Perturbation Theory (HH%PT) 15Hl7j]. in which the charmonia are treated non 



relativistically and coupled to the D mesons and Goldstone bosons in a manner consistent 
with heavy quark and chiral symmetries. In a non-relativistic theory the D meson kinetic 
energy scales as m^v 2 and momentum scales as mov, where mo is a D meson mass and 
v ft 1/2 is the typical velocity of the D mesons in the loops. With this scaling, Ref. [yj] 
showed that the loop diagrams with D mesons should be enhanced over tree level couplings 



1 The decay rate ip' — > rj c j vanishes in the non-relativistic quark model due to the vanishing overlap of the 
orbital wavcfunctions of the ip' and the ?y c , and is no longer zero once relativistic corrections are taken 
into account. However, quark_ models that include relativistic corrections still have trouble reproducing 
the correct rate for ip' — > 7/7 
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by a factor of 1/v. The rates for ip' — > J/ipn and ip' — > J/ipt] are sensitive to the product 
(72(725 where the J /if) coupling to D mesons is #2 and the ip' coupling to D mesons is g 2 . 
Ref. [l4| found a value of g 2 g' 2 consistent within errors with the two experimentally measured 
rates. This resolves the disagreement between the value of m^/m^ extracted from these 
decays and other extractions, since the prediction for the ratio of rates in terms of m u /m^ 
relied on the rates being dominated by the tree level HHxPT coupling. The value of g2g' 2 
extracted by Ref. [l^] is consistent with power counting estimates of g^ and g' 2 , which are 
both expected to be ~ (m c t> c )~ 3 / 2 up to constants of order unity. Other hadronic and 
radiative charmonium decays are also analyzed within the same formalism in Refs. 18N20|. 

The goal of this paper is to apply the same theory to the radiative decays mentioned 
above. One of our aims is to check whether the theory can also successfully resolve puzzles 
in radiative decays as one would hope. It is also important to check that couplings extracted 
from the hadronic decays are consistent with data on radiative decays. An important aspect 
of our analysis is that unlike Refs. Q, [l^, tree level counterterms are included in our 
calculations of both hadronic and radiative decays. Ref. 18] argued for an additional factor 
in the loop graphs of 1/(Attv j?) ~ 0.5 — 0.6, for v 2 C ~ 0.25 — 0.3, which would compensate the 
1/v enhancement of the loops. This factor, and the fact that v is not very small, support 
including both the loops and tree level interactions in the calculation, which we will do in 
this paper. This can have an important impact of the extracted values of the couplings #2 
and g 2 . Finally, an additional motivation for our analysis is that the extracted couplings are 
important for the physics of the X(3872) and other recently discovered charmonium bound 
states that have been interpreted as charmed meson molecules. If the X(3872) is a charmed 
meson bound state, then the coupling g 2 (g' 2 ) is an important theoretical input for calculations 
of X(3872) — > J j +X, so extraction of g 2 and g 2 is relevant to unconventional as well as 



conventional charmonia. For theoretica 
using effective field theory, see Refs. 121 



calculations of X(3872) to conventional charmonia 
Ll-23]. 

Our main result is that in order to obtain a consistent fit to both radiative decays as 
well as the hadronic decays considered in Refs. [jj], counterterm contributions must 
be included and the values of g 2 and g 2 will then be smaller than estimated in an analysis 
containing only the loop diagrams by a factor of 2.3. This decreases the overall size of 
the loop amplitude by a factor of 5. It is not possible to get reasonable agreement with 
radiative decay data without including counterterms. Since NRQCD is the microscopic 
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theory of charmonia, and does not include loop effects from charmonia, one is tempted to 
identify the result of a calculation of the J/ip — > rjc'-f amplitude in NRQCD with the tree 
level coupling in HH^PT. This is somewhat tenuous as the bare coupling in our theory 
has an infinite piece that must cancel the linear divergence in the meson loop integrals. 2 
Nevertheless, we regard it as satisfying that the size of the counterterms we extract in our 
fit with the smaller values of #2 and g' 2 are consistent within a factor of 2 with the quark 
model and NRQCD calculations of the radiative transitions. For other extractions of the 



couplings #2 and g' 2 in different theoretical frameworks, see, e.g., Refs. 24j-|30(. In Refs. 



31 



32] . the charmed meson loop corrections to radiative J/ip and ip' decays are studied in 
a version of HHxPT with relativistic propagators and couplings, as well as form factors at 
the vertices that regulate ultraviolet divergences. The form factors introduce an additional 
parameter into the calculations. These authors did not attempt to simultaneously fit the 
hadronic decays but used values of g% and g' 2 consistent with those obtained in Refs. 
Their results are also consistent with the experimental data on the radiative decays 
The effective HH^PT Lagrangian relevant to the hadronic decays is 



21 



33 



C = Tr[Hl (id + H a ] + ^Tr[HloH a a} - 9 -Tr[HlH b a ■ u ab ] 

+ ij (Tr[JV Jt] - Tr[,Fa\J'}) + i^Tr^HJ • V ' H a ] + H.c. . 



(1) 



Here H a = V a ■ a + P a and H a = —V a -cr + P a are the charmed and anti-charmed meson 
multiplets with V a and P a denoting the vector and pseudoscalar charmed mesons, respec- 
tively, and = if^ ■ a + r]c denotes the charmonium multiplets with -0^ anc i ^\ The a 
are the Pauli matrices, a and b denote flavor indices, and A d B = A(dB) — (dA)B. The 
first two terms in Eq. ([1]) are kinetic terms for the charmed mesons, A = mp* — is the 
hyperfine splitting, and m£>(m£)*) is the mass of pseudoscalar (vector) charmed meson. The 
third term contains the interactions of D mesons with the Goldstone boson fields which are 
contained in u = exp(i(f) / V2F) where <fi is a 3 x 3 matrix of Goldstone boson fields and 
F = 92.4 MeV is the pion decay constant. There are identical terms for the D mesons 
which are not explicitly shown. The tree level couplings for ip' — > J/ipir°(r]) come from the 
term with coupling constant A. The factor \- is defined by \- = u^x u ^ ~ U X^ U > where 
X = 2B ■ diag(m u ,m d ,m s ), m u m d and m s are the light quark masses and B = |(0|gg|0)|. 



This linear divergence is absent in dimensional regularization. 
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Finally, the IS charmonia couple to the D mesons via the last term with coupling g%. The 
same term, with J and g^ replaced with J 1 and g' 2 , couples the 2S charmonia to charmed 
mesons. 

The tree level decay amplitudes are 181 ]. 



iM{ip' -»■ J/^vr°) 



iAAe ijk qie^' ej 4 ' B du 
i(S/V?>)Ae l3k q4e J k ^B t 



(2) 



where 5 du = ^(m d 



m,. 



and 5,; = 4?(m 



m u +m d ■ 
2 



To leading order in the chiral 



1^0 



>sl — F 

expansion, these factors may be expressed in terms of light meson masses: B^ u - 
m 2 K+ + m 2 K+ — m n o)/F and B s i = (3/4) (m 2 — m^ o )/F. The 7r° — 77 mixing must also be 
included, and the mixing angle is 



1 ra^o - m 2 K+ + 



x/3 — ^ 

V 77 7TQ 



(3) 



When this mixing is included the first matrix element in Eq (T5]) is multiplied by 3/2. The 



18] . Since the 



loop diagrams contributing to the decay have been evaluated in Refs. 
decay to 7r°(?7) vanishes in the isospin (£77(3)) limit, the diagrams cancel in the sum over 
D°,D + , and appearing in the loop in the limit that all these mesons are degenerate. 
Mass differences between the mesons render the cancellation incomplete and are responsible 
for the finite contribution. 

For electromagnetic decays, we need to add couplings to the magnetic field and gauge 
the interactions in Eq. ([1]) . The tree level coupling of the charmonia to the magnetic fields 



is given by 34 



P 



^Tr[JB ■ a,P] + P -(Tr\.J' B ■ a.P] + H.c.) + ^Tr[J'B ■ a. J 1 



l> 



(4) 



where B is the magnetic field. Due to the presence of Pauli matrices these terms break heavy 
quark spin symmetry. The first term is responsible for the decay J ftp — > rjcj, the second for 
"0' — >* Icli and the third for ip' — > rfj. For the loop corrections to the radiative decays, we 



must a 



by [33 



so include the coupling of the charmed mesons to the magnetic field, which is given 



36] 



^-Tr[HlH b a ■ BQ ab ] 



2m, 



■Q'Tr[Hl& ■ BH a ] 



(5) 



where Q a b = diag(2/3, —1/3, —1/3), Q' = 2/3, and m c is the mass of charm quark. These 
terms are responsible for the decays D* — > Dj. Including leading as well as AQ C D/ m c 
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FIG. 1: Triangle diagrams with intermediate charmed meson loops. The charmed meson couplings 
to photon come from Eq. ([5]). 
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FIG. 2: Triangle diagrams with intermediate charmed meson loops. The charmed meson couplings 
to photon come from gauging the kinetic terms in Eq. ([TJ). 



suppressed terms is crucial for reproducing observed D* — > rates 36J. Ref. 33[ finds 
that a good fit to the experimental rates is obtained for the values m c = 1.5 GeV and 
P = 3.0 GeV" 1 . 

These couplings enter the radiative decays of charmonia through the triangle diagrams 
shown in Fig. [TJ There are also interactions that arise from gauging the derivatives in 
Eq. ([TJ. Gauging the derivatives in the kinetic term for the D mesons leads to couplings 
to the photon which contribute to the radiative decays via triangle loop diagrams shown 
in Fig. [2j Gauging the coupling leads to a contact interaction that directly couples 
charmonia, heavy mesons and the photon field, which is given by 



eg 2 Tr[J^H a a-AH a ]+H.c. 



(6) 
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where a = 2 or 3 only, i.e., only charged and strange D mesons appear in the interaction 
term of Eq. ([6]). The loop diagrams with contact interactions are shown in Fig. |3j 
The tree level amplitude for the J/ip — >■ rj c ^ decay, for example, is 

iM = pe ijk q i e]e J k '' 4 ' , (7) 

where q denotes the momentum of photon, and e 7 and e J ^ are the polarization vectors of 
the photon and J /if), respectively. The corresponding decay rate is 



71 m j/i, 

= -^^\qf. (8) 

Here the factor ^JfnJJ^m^ comes from the normalization of nonrelatvistic fields in HHxPT. 
In the non-relativistic quark model, p = 2ee c /m c « 0.3 GeV -1 , and the predicted decay rate 
is about a factor of 2 too large. The contribution to the amplitude from meson loops is also 
proportional e i j k q i e l -e'jj'^ so the loops give an additive shift to the p term for each decay. We 
evaluate the loops in pure dimensional regularization so linear divergences do not appear 
and the corrections from all loops are finite. The full rate is Eq. (JS} with Aio replaced by 
•Mfuii, where the Aifuii includes both the tree level interaction and the contributions from 
neutral, charged, and strange meson loops. The explicit expression for M. f u ii can be found 
in the Appendix. 

Before proceeding to our fits to the data, we will briefly discuss the power counting for 
the diagrams we have shown. As stated earlier, for non-relativistic D mesons one takes 



E ~ m^v 2 , p ~ ttidv, so the propagators scale as (mct) 2 )^ 1 and the loop integration 
measure is mj^v 5 . The vertices coupling the charmonium to D mesons carry a factor of 
p ~ m£>v . To estimate v 2 one may take the difference between an external and two internal 
mesons, so m D v 2 = |m charmonium - m Dpaix \ where m charmonium is the mass of one of the 
external charmonia and m£) pair is the mass of two D mesons in the loop. This leads to 
an estimate ranging from v 2 = 0.09 (m c harmonium = tti^ and m^pair = m-D + rn D *o) to 
v 2 = 0.5 (^harmonium = mj/y and mD pa ,ir = 2m D *o). Naively with this counting the triangle 
diagrams scale as (mj :) v 5 )(mDV 2 )~ 3 ('m£)v) 2 q = m 3 D vq. The first factor comes from the loop 
integration factor, the second from the propagators, the third factor from the derivative 
couplings of charmonium to D mesons, and the factor of q is the photon or pion momentum 
which comes from the coupling of these particles to D mesons. The diagrams with the 
contact interaction scale as (m%v 5 )(mov 2 )~ 2 q = m 3 D vq which is the same as the triangle 
graph. This is because there is no derivative in the contact interaction with the photon, the 
derivative in the charmonium D meson coupling must turn into the factor of q required by 
gauge invariance or chiral symmetry and there are only two propagators. The factor q is 
common to all diagrams including the tree level diagrams. Factors of mo are compensated 
by other dimensionful couplings so we will focus only on counting powers of v from here on. 
So the triangle graphs and graphs with the contact interaction are v suppressed relative to 
the tree level interactions. However, these hadronic decays violate either isospin or SU(3) 
and the radiative decays violate heavy quark symmetry so there are cancellations between 
graphs due to heavy meson mass differences that are missed by this power counting. One can 
formally modify the power counting counting in the following way. The inverse propagator 

2 

for a non-relativistic meson can be written as E — tt — h b + 5, where the residual mass term 
in the propagator has been split into a term b which is common to all D meson states and 
a term 5 contains SU(3) breaking and hyperfine splittings that are different for different D 
mesons. Expand the D meson propagators as 

1 _ 1 5 

£-& + ^ = £-^~( £ -^ + 6) 2+ "" 

The graph in which all propagators contribute only the first term is zero by symmetry. 
In order to get a non-vanishing result at least one propagator in the graph must give a 
contribution from the second term, then the power counting says the graph is enhanced by 
a factor of S/m D v 2 , which makes the graph 1/v enhanced rather than v suppressed relative 
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to the tree level diagrams [jj], [lsj]. Since v is not very small this could be compensated 
by other numerical factors. In practice it is easier to simply calculate the graphs with the 
unexpanded propagators but expanding the propagator makes it clear that after summing 
over all graphs one gets a 1/v enhancement. In this paper, we will take the viewpoint that 
the leading one loop diagrams are of roughly the same size as the tree level contributions 
and include both in the decays, then try to simultaneously fit the radiative and hadronic 
decays mentioned above. 

A separate question is whether higher order chiral corrections are under control. Certainly 
some chiral corrections are suppressed as argued for different charmonium radiative decays 
in Ref. {20)]. But in a subgraph with a ladders of single pion exchanges between a pair of 
D mesons, non-relativistic power counting shows that the ladder with n + 1 single pion 
exchanges is suppressed relative to one with n pion exchanges by a factor g 1 mr>vl {8ttF 2 ) = 
pj (320MeV) 37J where p is the relative momentum of the D mesons. This would require 
p = rri£>v with v < 0.08 to be less than 1/2. In some channels a resummation of single pion 
exchanges may be needed to do accurate calculations. Such a resummation is beyond the 
scope of this paper. Here we are simply interested in the impact that including the tree level 
interactions and simultaneously fitting the radiative and hadronic decays has on the values 
of g 2 and g' 2 and therefore the size of D meson loop contributions to charmonium decays. 

To constrain the parameters g 2 and g 2 , we determine the parameter A and the product 
g 2 g 2 from the measured rates for T[ip' — > J/tpn ] and Tty' — > J/iprj\. Because the predictions 
for the decay rates are quadratic in g 2 g' 2 , this does not completely determine g 2 g' 2 , but yields 
two possible solutions. Then we fix the relative size of the two couplings using the relation 



92 — 92y m 'ip/' m J/^i which follows if the dimensionless coupling of the J/ip and ip' to D 
mesons is the same Once g 2 and g 2 are determined this way from the hadronic 

decays, the only parameters remaining in the radiative decays are p, p', and p", which can 
be determined from the three decay rates Y[J /ip — > rjcj], Tfy' — > ^7], and T[ip' — > rj' c j]. The 
results of determining A, g 2 , and g' 2 are shown in the first three columns of Table 1. One 
possible fit to the hadronic decays yields A = — 6.36 10~ 3 GeV -2 and g 2 g 2 = 1.73 GeV~ 3 . 
This is a very small value of A, almost two orders of magnitude smaller than the estimate 
A ~ l/(2m 2 c ) in Ref. 18]. This fit yields a value of g 2 g' 2 similar to that of Refs. Q, \\\ . 
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In this case, we get a value of 17252 that is a factor of two smaller than Refs. because our calculati 



10ns 



9 



A {GeV' 2 ) 


g 2 (GeV~ 3 / 2 ) 


g' 2 (GeV~ 3 / 2 ) 


p (GeT/- 1 ) 


p" (GeF- 1 ) 


p' (GeV- 1 ) 


-0.00636 +4 ^ 

— 4X1U 


-i 07+0.02 

J.-'J ' —0.02 


1.26i°;^ 


1-88185S 


0.83±°il 


1 5fi+ a05 

—0.05 


-0.00636+JSJi; 

— 4X 1 U 


1 07+O.O2 
±•01 _q.02 


1 26 +a02 
—0.02 


9 97+O.O9 
—0.09 


i.26±g:« 


1 50+0.05 


u - uzo ' -0.0007 


n ^qq+0.023 
u.djy_ 023 


u.o*y_ 021 


o-i9i±g5iS 


-0.0192^ : ^| 


n OO7+0.024 
u ' zo ' -0.024 


u - uzo ' -0.0007 


n ^qq+0.023 
u.djy_ 023 


n 540+0.021 


0.598+°;^ 


n41 r+o.i05 

U.41O_ 1Q5 


n 01 q+0. 024 
u - o±o -0.024 








0.275 


0.263 


0.0417 



TABLE I: The numerical results for fitting parameters to hadronic and radiative charmonium 



decays. The fit is explained in the text. The quark model 



9j] predictions for the parameters p, p' 



and p" are shown in the bottom line. Errors are due to experimental uncertainties only. 



The second possible fit is A = 2.57 10 -2 GeV -2 and g 2 g 2 = 0.329 GeV~ 3 , a value 5.3 times 
smaller than the first fit. The value of A is closer to the estimate of l/(2m 2 ), but still a 
factor of 10 smaller. The results of fitting the parameters p, p', and p", are shown in the last 
three columns of Table 1. For each choice of A, g2, and g' 2 , there are two possible solutions 
for p, p', or p", for a total of four possible solutions. The values of p and p" are much closer 
to the quark model predictions (shown in the bottom row of Table 1) in the fit with a smaller 
value of (?2fl , 2- The extracted value of p' does not come close to the quark model prediction of 
Ref. {9 1, but this model does not give a good prediction for the rate ip' — > r\ c ^f. In the fits with 
the larger value of gig 2 the extracted values of p, p', and p" are much larger. This indicates 
that for these choices of parameters fine tuned cancellations between the tree level and loop 
diagrams are required to fit the data. This can also be clearly seen in Table II, where we 
give the loop contribution to the decay for each fit. For the first fit with g 2 g' 2 = 1.73 GeV -3 , 
T[J/ij} — > rjcj] is over predicted by a factor of 100 and T[ip' — > 77^.7] is over predicted by a 
factor of 20 without the counterterm contribution. Thus, in order to fit these decays, fine- 
tuned cancellations between loop and tree level contributions must occur. Though the loop 
contributions by themselves do not do a good job of producing the radiative decay rates for 
the smaller value of g2g^ the discrepancy is not nearly as large. 



of the loop amplitudes for ip' — > J '/ 't/jn° '(77) disagree with the analytic results of Refs. [3, H] by an overall 
factor of two. This is because we include graphs in which the 7r° or r\ couples to the &•*> mesons, instead 
of the £>W 

mesons, that arc omitted in Refs. 

HQ. 
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Fit 1 


Fit 2 


QM 


PDG, BES III, and CLEO 


r[J/V> ->■ ??c7];oo P 


163^° keV 


5.96lg;|J keV 


2.9 keV 


1.58± 0.37 keV[10] 


-»• r]' c % oop 


3.30±°;|S keV 


0.119±g : °l? keV 


0.21 keV 


0.143± 0.027±0.092 keV[38] 


T[ip' -)• r/ c 7] ioop 


16.34±g;|| MeV 


597±|? keV 


9.7 keV 


0.97± 0.14 keV[39j 



TABLE II: The decay rates contribution from the loops alone for the two solutions for g 2 g 2 * s 
compared with the results in the quark model [9(] and experimental data. Fit 1 corresponds to 

9292 



- 1.73±g:g| (GeV" 3 ), Fit 2 corresponds to g 2 g' 2 = 0.329^;^ (GeW 3 ) 



For the decay if)' — > ri c 7 there are rather severely fine tuned cancellations between loop 
diagrams and tree level contributions for both fits. The photon energies in the decays 
J /if) — > rjc^, if)' — > ?7 C 7, and if)' — > 7^.7 are 114 MeV, 638 MeV, and 49 MeV, respectively. The 
photon energy in the second decay may be too large for either the quark model or low energy 
effective theory to be accurate. As an alternative approach, one can simply try to extract g 2 
and g' 2 independently from the radiative decays J /if) — > 7^7 and if)' — > 77^7, using the quark 
model js] to estimate the parameters p and p" . Since the decay J /if) — >■ r/ c 7 is quadratic 
in g 2 and the decay if)' — > 77^.7 is quadratic in g' 2 2 , there are two possible solutions for each 
parameter. We find g 2 = 0.255±g;g^ GeV ~ 3 / 2 or 0.659±g;gi| GeV" 3 / 2 and g' 2 = 0.264t°; 2 ^ 
GeV" 3//2 or 0.855lQ;og[} GeV" 3//2 . Note that the value of g 2 g 2 obtained this way is also smaller 
than the value obtained in the first fit to the combined hadronic and radiative decays. Also 
the ratio g 2 /g 2 obtained using the the smaller two central values is 1.04 while using the 
larger two central values the ratio is 1.3. Both of these are a little larger than one expects 



from the hypothesis g' 2 /g 2 = ^/rnj/^/m^' = 0.92. 

In summary, we have computed the decay rates for J /if) — > r] c ^ and if)' — > r/^'y including 
tree level and one loop diagrams with charmed mesons in HHyPT. We combined our results 
with the decay rates for if)' — > J/if)Tx Q {j]) found in Refs. 14], isj], used the relationship 
g' 2 = g 2 y/mj/^/m^i, and fit the five remaining coupling constants simultaneously. Including 
tree level couplings is essential for simultaneously reproducing all the decay rates. A smaller 
value of g 2 g' 2 = 0.33 GeV -3 is required to avoid large cancellations between tree level and 
charmed meson loop contributions to the radiative decay. The tree level couplings p and 
p" in this fit are consistent (to within a factor of 2) with expectations based on the quark 
model and NRQCD. 
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APPENDIX 

In this Appendix, we calculate the loop diagrams for the Ml radiative decays. We 
present the calculation of J/ip — > rj c ^, the calculations of the decays ip' — > rj^ (-0' — > 7^,7) 
are obtained by replacing rrij/^ with m^> and g\ with g 2 g' 2 (g? 2 )- 

The triangle loop diagrams in Fig. [1] have a similar form as the triangle loop diagrams in 
the hadronic decays, therefore the notation used here will be the almost the same as that 



of Ref . 18] . We refer the reader to that paper for explicit expressions for the integrals. The 
amplitude from Fig. [TJ^a) is 

_dH e ijkgi e](2l - q) k e J ^ ■ I 

(2 " )4 8(*o - ^ + ie)(l + ^ + b DD -ie)(l -q -^-A + re) 



iMia = -2igj\ m J 



where boD = 2bid — mj/^, Ai = |(e/3 + is relevant for loops with neutral D mesons 
and A3 = — |(e/3 — ^-) is relevant for loops with charged and strange D mesons. Here 
l[ 2 \q, mi, m 2 , m 3 ) only differs from the function defined in Ref. [3] by omitting a factor of 
mim 2 m^ from the denominator. Fig. QJb) contributes 



-j e k \Q\ 

x (2lf\q,m D ,m D *,m D ,) + Al[ 2) (q, m D , m D *, m D *) - I {l) (q,m D ,m D ,,m D *)) 



iMu = 2g 2 X m e ijk q i e]e k / \q\ (10) 
r(2)/„ ™ ™ m \ , atV), 



where the functions Ig 2 \q, mi,m2, m^) and h(q, mi, m 2 , 7713) are again the same as functions 



in Ref. 



up to a factor of mim 2 m 3 . Fig. []Jc) contributes 



iMi c = 2g 2 2 X m e tjkqi e]e J k ^\q\ 2 (11) 
x (2I^\q,m D ,,m D *,m D ) + Qlf\q,m D *,m D »,m D ) - I {l \q,m D *,m D ,,m D )) . 
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Fig. QJd) contributes 

iM ld = 2g 2 2 \ m e ljk q i e]e J k / " p \q\ 2 (12) 
x (2I^\q,m D *,m D ,m D *) + 4l[ 2 \q,m D *,m D ,m D *) - I w (q,m D »,m Dl m D *)) , 

where A2 = — |(e/3 — ^-) is relevant for loops with neutral D mesons and A4 = |(e/3 + ^-) 
is relevant for loops with charged and strange D mesons. Finally, Fig. [1(e) gives 



t j t k \q\ 

x (2I i 2) (q,m D , J m D ,,m D ,) + 8l[ 2) (q, m D * : m D *, m D ,) - I {1 \q,m D *,m D ^m D ,)) 



iMi e = 2g 2 X 2{A) e ijk q i e]e k / \q\ (13) 
r( 2 )f„ ™_ ™_ ™_ Ufiff 2 )/ 



In addition to the triangle diagrams with the couplings of D and D* mesons to the magnetic 
field from Eq. 03), there are also two triangle diagrams with the coupling of the photon to 
charged D and D* mesons that arises due to gauging their kinetic terms. These are shown 
in Fig. [2j The sum of these two diagrams yields 

iM 2 = "igWijkqi^^Wl 2 ( — Ii\m D ,rn D *,m D ) — l[ 2 \m D *, m D , m D *) ] . (14) 

J \m D m D * J 

So far we have only included the interactions coupling the photon to D and D* mesons. 
There are additional diagrams where the photons couple to D and D* mesons that give an 
equal contribution. 

Fig. [3] shows the loop diagrams with the contact interaction that arises from gauging the 
coupling g 2 . Fig. Eta) yields 



iM 3a = i2g 2 e J 



dH e^l + q^eje^ 



(27F)4 4(^o - i + ie)(lo + So + §^ + b DD . - re) 
-glee i j k q i e! j E k Jli 'T '(m D , m D *) , (15) 



where 



l'(m D ,m D *) = ( mD mD ) J — \q\ 2 + 2/i DD *(b DD * + q ) , (16) 

An \m D * + m D J y m D + m D * 

for fiDD* = m D m D * / '{m D + m D *) and boD* = m D + m D* — m J/ii>- Fig. Mb) is related to 
Fig. G^a) by charge conjugation so 

iM 3b = iMsa = -glee ijk qie]e J k /n 'I f (m D ,m D *) . (17) 
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The graphs in Fig. [3]^c) and Fig. |3]^d) both vanish, so the total contribution to the amplitude 
from loops with contact interactions is 2iM.- ia . Only diagrams with charged and strange D 
mesons in the loop will contribute. 

The total amplitude from loop diagrams in Fig. [T] with neutral D mesons is given by 

iM\ = e ijk q}e]e^{gl\q) 2 [4\il[ 2 \q,m D ,m D ,rn D *) + 2\ l (2I {2 \q, m D , m D *, m D *) 

+4I < L 2 \q,m D ,m D *,m D *) - I w (q,m D ,m D *,m D *)) + 2\ 1 (2I^\q,m D *,m D *,m D ) 
+QI < L 2 \q,m D *,m D * 7 m D ) - I {l \q,m D ,,m D *,m D )) + 2\ 2 {2I^\q,m D *,m D ,m D *) 
+4l[ 2) (q,m D *,m Dl m D *) - \q,m D * ,m D ,m D *)) + 2\ 2 (2l{ ) 2 \q,m D * 1 m D *,m D *) 
+8l[ 2 \q,m D *,m D *,m D *) - I {1) (q,m D *,m D *,m D *))]} . (18) 

Adding the five diagrams with the photon coupling to a D or D* doubles this contribution. 
The contribution from diagrams of Fig. [1] with charged and strange charmed mesons in the 
loops is obtained by substituting Ai with A3 and A2 with A4. In addition, the contributions 
from the diagrams in Fig. [2] and Fig. [3] with charged and strange charmed mesons in the 
loops need to be included. The full decay amplitude is 

iMfuii = iM Q + 2{iM\ + %M\ + %M\ + iM c 2 + %M\ + %M% a + iM s 3a ) , (19) 

where the superscript c(s) indicates a contribution from loops with charged (strange) D 
mesons. The decay rate for J/?/> — > jrj c is given by 

T[J/^ -)■ 1Vc ] = ^(V^W%M /(lll ) 2 ^. (20) 
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